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Abstract. 

We derive the master equation of a system of two coupled qubits by taking into 
account their interaction with two independent bosonic baths. Important features 
of the dynamics are brought to light, such as the structure of the stationary state 
at general temperatures and the behaviour of the entanglement at zero temperature, 
showing the phenomena of sudden death and sudden birth as well as the presence of 
stationary entanglement for long times. The model here presented is quite versatile 
and can be of interest in the study of both Josephson junction architectures and cavity- 
QED. 
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1. Introduction 

During the last two decades the problem of controlling the coherent time evolution of 
quantum systems has received a great deal of attention both because it could lead to 
a deeper understanding of fundamental aspects of quantum mechanics and because 
it is of interest for applications pQ. The possibility of generating and controlling 
non classical correlations (i.e. entangled states) in multipartite systems, despite their 
coupling with an external environment, is a fundamental ingredient for instance in the 
theory of quantum measurement and in the study of the border between the quantum 
world and the macroscopic classical world [2]. In addition such a possibility is an 
essential goal in the field of quantum computing and quantum information theory [3]. 
It is well known that the interaction of an open quantum system with an external 
reservoir is an important source of dissipation and decoherence. In order to describe such 
phenomena, a master equation approach can be used [HOE]. In particular, following 
Ref. [6], one has at his disposal a general formalism allowing to derive the quantum 
master equation of a general quantum system, provided one knows the Hamiltonian 
governing the unitary part of its dynamics. Exploiting this recipe, one can see that 
the dissipative dynamics, in the Born-Markov and rotating wave approximations, is 
described by a master equation in which the quantum jumps occur among eigenstates 
of the Hamiltonian of the open quantum system under scrutiny. This microscopic 
approach may be sometimes in contrast with some phenomenological approaches to 
quantum dissipative dynamics present in the literature, as one can see for example in 

Refs. mm!. 

Within this framework and with the aforementioned approach, here we analyze the 
dynamical behavior of the entanglement between two coupled two-state systems each of 
them interacting with a bosonic bath. The model is quite versatile and can be exploited, 
for instance, in order to describe the dipole-dipole (flux or charge) interaction of two 
distant atomic (flux or charge) qubits. In the case of Josephson flux qubits the coupling 
term corresponds to a flux-flux coupling proportional to their mutual inductance [10]. It 
is worth noting that the same model can be used to describe richer physical situations, 
for instance the coupling between a Josephson junction based qubit in the charge regime 
with an impurity located in the substrate of the device [13] . In this case the qubit can 
be thought as coupled with external environmental degrees of freedom describing the 
auxiliary circuitry (necessary, for instance, for the qubit control and the redout) which 
can be modeled as an infinite bath of harmonic oscillators, while the substrate impurity 
may interact with a phononic bath. 

Starting from such a model of qubit-qubit interaction, and including counter- 
rotating terms in the system Hamiltonian, we derive a quantum master equation in 
order to describe the dissipative dynamics of the two-qubit system, concentrating on 
some aspects such as the structure of the stationary state of the system and the time 
evolution of the entanglement between the two qubits. The entanglement dynamics of 
open bipartite quantum systems has been studied in previous works, for example in Refs. 
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[TTJ [T2J [33], bringing to light important features such as the complete disentanglement 
of the system in a finite time (entanglement sudden death), and, in the case of Ref. 
[H], the sudden reappearance of entanglement (sudden birth). In the latter case, the 
presence of stationary entanglement for very long times has been related to the presence 
of a reservoir common to the two subsystems, so that the non-local correlations between 
them can be thought of as due to the mediation of the environment. In this paper we 
will show that stationary entanglement can occur also in the case wherein the two 
qubits interact with two independent reservoirs, provided the interaction between the 
two subsystems contains also counter-rotating terms, which are usually neglected in 
the literature. In such a case, the stationary entanglement occurs because the counter- 
rotating terms cause the ground state of the system to be an entangled state of the 
bipartite system. 

The paper is structured as follows. In section 2 we introduce the model previously 
described and we microscopically derive the Markovian master equation in the weak 
damping limit assuming that the two reservoir are independent and with arbitrary 
temperatures, 7\ and T 2 . Therefore, in section 4, we analyze the dynamics of the 
system in the limit case T\ = T 2 = discussing the stationary entanglement and the 
features of its time evolution considering different initial conditions for the bipartite 
system. Finally, conclusive remarks are given in section 5. 

2. The model 

Let us consider two interacting two- level systems and let us call |0)i ((0)2) the ground 
state of the first (second) system and |l)i ( 1 1) 2 ) the corresponding excited state. Let us 
assume that the two systems are coupled in such a way that their unitary dynamics is 
governed by the following Hamiltonian (in units of K): 



where uji is the Bohr frequency of the z-th two- level system, A/2 is the coupling constant 
and where we have used the Pauli operators 0+ = |1}^(0|, <r_ = |0}jj(l| and 
o~x = 0+ + a_, with % = 1,2. In the case of Josephson flux qubits the coupling 
term in Eq. ([TJ) corresponds to a flux- flux coupling with A/2 proportional to their 
mutual inductance [TO] . 

It is worth noting that in the Hamiltonian also the counter-rotating terms of the 
interaction have been included and we will see that they play a central role in the 
dynamics of the entanglement between the two systems. 

The model in Eq. (CQ) can be exactly diagonalized. By exploiting the fact that 
the Hamiltonian Hs in the uncoupled basis {|00), |11), |10), |01)}, where for instance 
1 00) = 1 0) 1 <g> |0)2, is block diagonal, it is straightforward to show that the eigenvalues 
(given for increasing energies) are: 



H S = w®*® + ^1V 2) + ± aVa®, 



(1) 
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E b = - {uji + uj 2 ) - ^ \J (^2 - ^i) 2 + A 2 



E c=\ (wi + ^2) + ^yVa-c^f + A 2 



£ rf = -K + cu 2 ) + (^ 2 + ^)" + A 2 , (2) 
while the corresponding eigenstates are: 
I a) = cos — 100) - sin— 111) 

1 6) = cos— 1 10) - sin— 101) 

|c) = sin 1 10) +cos^|01) 

\d) = shAoO) + cos^|ll). (3) 
Here oj 2 > ^1 and the parameters and #77 satisfy the relations: 

sint// = — , cost'/ = — (4 J 

y (cj 2 + cui) 2 + A 2 y (cj 2 + ^i) 2 + A 2 

sm0 n — — , cos On — — j^^^=^^^= (5) 

y (uj 2 - ^i) 2 + A 2 y (cj 2 - cji) 2 + A 2 

The losses in the system under scrutiny will be taken into account by considering 
the coupling between the «-th system and its own reservoir at temperature Tj. In 
the following we will consider the case of independent bosonic reservoirs, whose 
temperatures, in the general case, can take different values. The total Hamiltonian 
of the bipartite system and the reservoirs can thus be written as follows: 

H — H s + H E + H int 



k j 

H illt = o-^ ®J2 tk ( ak + a l) + a * ] ® J2 9i ( b J + & 1) • ( 6 ) 

k j 

Here (a' k ) is the annihilation (creation) operator of the fc-th mode (of frequency 
Uk) of the reservoir interacting with the first subsystem and similarly bj (bj) is the 
annihilation (creation) operator of the j-th mode (of frequency uOj) of the reservoir 
interacting with the second subsystem. The parameters and gj are the corresponding 
coupling constants. 
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From this model it is possible to microscopically derive the Markovian master 
equation describing all the relaxation phenomena in the dynamics of the bipartite system 
under study. To this end in the next section we will exploit the general formalism given 
in Ref. [6j. 

3. Derivation of the Markovian master equation in the weak damping limit 

From the model in Eq. ([6]) it is possible to microscopically derive the master equation 
for the evolution of the bipartite system, by following the general procedure outlined in 
Ref. [6J , in the Born-Markov and rotating wave approximations^. The main point of the 
formalism is that all the jump processes involve transitions between dressed states of 
the open system under study, i.e. the eigenstates of the Hamiltonian H s of the system, 
which in our case is given by Eq. (jHJ). 

It is well known that, in the Schrodinger picture, the Markovian master equation 
for a generic open quantum system with Hamiltonian H is given by [6j: 

p{t)= -i[H ,p{t)\ 

+ EE^M U(w)p(t)4H - \ {4M4j("M*)}) » (7) 

ui a,/3 ^ ' 

where the symbol {•, •} denotes the anticommutator between operators. In the derivation 
of Eq. (J7|) we have to assume that the interaction Hamiltonian between the system and 
the environment is of the form Hi = A a <E> B a , where A a = A^ a acts on the Hilbert 
space of the open system under scrutiny, while B a = acts on the Hilbert space of 
the environment. In particular, in the case of the model in Eq. ([6]) the sum consists 
of two terms, since we have assumed that each two-level system interacts with its own 
reservoir. In Eq. (0), a renormalization Hamiltonian has been neglected and the rates 
7a,/? (k*) are given by the Fourier transforms of the environment correlation functions, 
according to: 

/+oo 
dre*" (Bi(r)B^O)) . (8) 
-oo 

Concerning the jump operators A a (u), their number is given by the number of different 
Bohr frequencies relative to H and they can be calculated from the relation [6j: 

A a (u) = n(e)A,n(e'), (9) 

e' — e=ui 

where n(e) is the projector on the eigenspace of the open system relative to the energy 
e and the sum is extended to all the couples of e and e' such that e' — e = uj. 

\ We stress the point that in this paper we will always call rotating wave approximation the operation 
of neglecting rapidly oscillating terms in the dissipative part of the master equation. In addition, 
we note that this should not be confused with the elimination of the counter-rotating terms in the 
Hamiltonian of the system which will always be taken into account. 
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Figure 1. Schematic illustration of the eigenstates of the bipartite qubit-qubit system 
and of the allowed transitions characterized by the Bohr frequencies u>i and loh. 



For our model we can identify the operators A\ = a x , A 2 = o~x , B\ = 
J2k e k ( a fc + a l) anc ^ ^2 — J2j 9j (bj + bjj . By calculating the matrix elements of A\ 
and A2 between the eigenstates of H$ given in Eq. (J3j) , it is possible to see that there 
are only two possible values for the Bohr frequencies of the transitions allowed. The 
first one is 



(U 2 + C0\) 2 + \ 2 - J (u 2 - UJlf + A 2 



for the transitions \b) — > \a) and \d) — > |c), corresponding to the jump operator: 

J n = (a|Ai|6)|a)(6| + (c|Ai|d)|c)<d| 
relative to the coupling of the first qubit with its own reservoir, and to: 

J I2 = (a\A 2 \b)\a)(b\ + (c\A 2 \d)\c)(d\ 

due to the coupling between the second qubit with the corresponding reservoir. 
Similarly, the second Bohr frequency is: 



(10) 
(11) 



]- (yj(u2 + ui) 2 + A 2 + yj (uj 2 - uJif + A ; 



L0 U 

for the transitions |c) — > \a) and \d) —>■ \b), corresponding to the jump operator: 

J in = (a\A 1 \c)\a)(c\ + (b\A 1 \d)\b)(d\ (12) 

relative to the coupling of the first qubit with its own reservoir, and to: 

Jin = (a\A 2 \c)\a)(c\ + (b\A 2 \d)\b)(d\ (13) 

due to the coupling between the second qubit with the corresponding reservoir. 

Inserting in Eq. ([7]) the structure of the jump operators given in Eqs. ( TTOj) - fTT3|) . 
the Markovian master equation can be cast in the following form: 

p{t)= -i[H s ,p{t)\ 
11 



/ 1 \ 

i=i 1=1 \ " ' 
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II 2 , 

i=i 1=1 ^ 



-fl.^)|] (in 



where, assuming that the two reservoirs are independent and each of them is in a thermal 
state, with temperatures T\ and T 2 respectively, one has: 

/+oo 
dre^/s/^S^O)), (15) 
-co 

and the Kubo-Martin-Schwinger relation 

Ti^e^^'Ti,,, (16) 

holds [6]. The hypotesis of independent reservoirs we have made consists in assuming 
that /s/(r)5 TO (0)\ = 0, i.e. 7 Mm = 0, when Z ^ m. 

Finally, by inserting Eqs. ( fT0l) -(fT3l) into Eq. ( fl4l) and rearranging the terms, we 
obtain: 

p{t) = -i[H s ,p{t)\ (17) 



+d \a)(b\p(t)\b)(a\ - - {\b)(b\,p(t)} )+ Cl j( \a)(c\p(t)\c)(a\ - - {\c) (c\, p(t)} 



+ CI [\b){d\p(t)\d)(b\ - - {\d){d\,p(t)}j + c H (Jc)<d|p(t)|d><c| - - {\d)(d\,p(t)} 
+cj (\b)(a\p(t)\a)(b\ -l{\a)(a\,p(t)}) + c H ( \c)(a\p(t)\a)(c\ - ]- {\a) (a\, p(t)} 



+cj {\d){b\p(t)\b)(d\ - - {\b){b\,p(t)}j + cu [\d)(c\p(t)\c){d\ - - {\c)(c\,p(t)} 
H-Ccr^ (|a><6lp(t)|d><c| H- lc><d[p(*)|6><a|) + Cc^j (Ia><c|p(t)|£l><6| H- |6><d|p(t)|c><a|) 



+c cr ,i (\d)(c\p(t)\a)(b\ + \b)(a\p(t)\c)(d\) + c crJI (\d)(b\p(t)\a)(c\ + \c)(a\p(t)\b)(d\) . 
The decay rates q, and the cross terms c cr ^ (with i — I, II), are given by: 

@i 9 ii , ■ ®i . ^ ii^ z 



- / = 7/,ii I cos — cos — + sin — sm — 

9 1 9j[ 9 j 9 n x 
+ li,22 [ cos — sm — + sm — cos — 

9 i . 9 u . 9 j 9 H ^ 2 
r , , = - , , j , \ cos — sin — - sm — cos — 

, 9 1 9 U . 9i . 9 n \ 2 
h ln,22 \ cos — cos — sm — sm — j , (19) 
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Ccr,I = 7/.11 I COS — COS — + Sill — Sill — 

/ 9 1 . On . . Oi H \ 2 / onN 
- 7z,22 I cos — sin — + sm — cos — I , (20) 

/ Oi . On . Oi On 
c cr ,ii = ~ 7//,n I cos -sm — - sm — cos — 

( Of 9 H . Oi . 0n\ 2 , 
+ lii,-22 I cos — cos — - sm — sm — I , (21) 

The corresponding excitation rates Cj, and the cross terms c cr>i) are obtained by 
substituting, in Eqs. (fT5j) - (l2"Tj) . j^u with the corresponding quantities 7i,i« : when the 
temperatures of the two reservoirs T\ and T 2 are both zero, all this coefficients vanish, 
according to Eq. (fT6]) . Physically this means that there is no possibility to create 
excitations in the bipartite system due to the interaction with the reservoirs. 

Let us conclude this section by reminding that the master equation has been 
derived in the framework of Born-Markov and rotating wave approximations. The last 
approximation is valid as long as the relaxation time of the system is much longer than 
the time characterizing its unitary dynamics [6]. Mathematically this is equivalent to 
assume that the coupling with the environment is weak enough to say that the relaxation 
rates are all much smaller than the smallest nonzero Bohr frequency relative to H$, that 
is uji. 

4. Dynamics 

The master equation given by Eq. (fT7|) . is equivalent to a system of coupled differential 
equation, the first four of which describe the time evolution of the populations of the 
dressed states \a), \b), |c) and \d), namely 



Paa{t) = - (Cj + Cn) Paa{t) + Cjp 66 (t) + C H p cc (t) 
Pbb(t) = Cip aa (t) - (Cf + C H ) Pbb(t) + Cnp dd (t) 
Pccit) = CllPaa(t) ~ (Cn + Ci) p cc (t) + Cip dd (t) 

Pdd{t) = c~nPbb{t) + cipcc(t) - (cj + cn) Pdd{t), (22) 

while the other equations describe the time evolution of the coherences: 

(cn + 2ci + cn) 



Pac(t) 



ILOn - 



Pac(t) + Cer,lPbd(t) 
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(23) 



(24) 



Pad{t) 



Pbc{t) 



1>V da - 



(Cj + C H + C! + C n ) 



lUJ cb - 



(c/ + c H + c/ + +c/j) 



Pad{t) 



Pbc{t), 



(25) 



where uJda = E d — E a and u) c b = E c — E^. The other coherences can be obtained by 
complex conjugation. In the following subsections we will discuss the existence of a 
stationary solution and we will bring to light the features characterizing the dynamics 
of the entanglement when both T\ and T2 are equal to zero. 

4-1. Solution I: The stationary state 

Let us look for the existence of the stationary solution by imposing pa = (Vz = 
{a,b,c,d}) in Eqs. fT22|) and the normalization condition Tr(p(t)) = 1. In such a 
condition, after some calculations, it is possible to prove the existence of a stationary 
solution given by: 

CiCu 



Paa, ST 



,ST 



Pec, ST 



{Cj + C/)(C// + C/j) 
CiCu 

(c 7 + C/)(C// + c n ) 
ClClI 

(C/ + C/)(C// + C/j) 
CiCu 



Pdd, st - 7 1 - w " — , - x (26) 

(c 7 + ci)(c n + c n ) 

As it is immediate to see, when T\ = T 2 = K, p aa ,ST — 1 ; namely the stationary 
solution coincides with the ground state of the bipartite system. In view of Eq. ([3]), 
this implies the existence of a stationary entanglement traceable back to the presence 
of counter-rotating terms in the interaction Hamiltonian (see Eq. ([T])) describing the 
coupling between the two two-state systems. We will discuss stationary entanglement 
in the next subsection. 



4-2. Solution II: Entanglement dynamics at zero temperature 



Let us consider now the analysis of the system dynamics when the temperatures of the 
two reservoirs, 7\ and T 2 , are equal to K. Moreover let us restrict our attention to 
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the case of exact resonance oj\ = 10% and in which the two environment have both flat 
spectra and are coupled with equal strength to the respective subsystems, which means 
that 7^11 = 7/7,11 = 77,22 = li, 22- From Eqs. (jSJ), (120"]) and ( |2T|) . it is straightforward 
to show that in this case the cross coefficients c cri 7 and c crj 77 are both equal to zero, so 
that each coherence between the dressed states of the system evolves independently for 
the other ones. 

Putting this condition into Eqs. (1221 - ff25l) . it is possible to show that their solution 

is: 

(27) 

vcn) \ (28) 
+Cll)t , (29) 

(30) 



and 



Paaif) 


= 1 - p bb (t) - p cc {t) 


-Pdd(t), 


Pbb(t) 


= (p bb (0) + p dd (0))e- 


- Clt ~ Pdd(0)e~ 


Pccit) 


= (pcc(0) + Pm(0)) e 


- Cnt -Pdd(0)e 


Pddit) 


= Pdd^)e- {CI+Cll)t 





Pacit) 
Pbdit) 
Pab{t) 
Pcd(t) 
Pad{t) 
Pbc{t) 



IUJJI- 



i c II+ 2c I+ c n) 



{ 2c i+ c ii+ c ii) 



20JJ-- 



( c I + 2c II+ c l) 



A, , ( C I+ C II + C I + C Il) 

LLU da 2 



. (ej + cjj + e/ + + c J - J -) 



PaM 

Pbd(0) 

Pab{0) 
Pcd(0) 

Pad{0) 
Pfcc(O) 



(31) 

(32) 
(33) 
(34) 
(35) 
(36) 



Starting from the initial state 



!#))) = V^|01) + e*y/T=j\10), 



p being a non- negative real number < 1, and by taking into account that 

/ f) \ f f) \ 2 

Pbb{0) 



1 — p cos 



On 
2 



^Jp sin 



On 
2 



Pcc(0) 
P6c(0) 



V 1 -p sin yfj + ^ cos (t") 



V 1 -p cos ( ~y) - Vp sin (~y 



x e 



-Up 



V 1 ~p sin [jfj + Vp cos 



(37) 

(38) 
(39) 



(40) 
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while pa = and pij = otherwise, it is possible to see that the density matrix describing 
the time evolution of the bipartite system at a generic instant of time t assumes the 
following form: 

p(t) = Paa(t)\a)(a\ + p bb (t)\b)(b\ + p cc (t)\c)(c\ + Pdd (t)\d)(d\ 

+ p bc (t)\b)(c\+p cb (t)\c)(b\ (41) 

To assess how much entanglement is stored in this bipartite quantum system at different 
instants of time we use the concurrence, a function introduced by Wootters [32], equal 
to 1 for maximally entangled states and zero for separable states, defined as: 



C(t) = max(0, y/\$j - - Vhit) ~ (42) 

where {Xi(i}} are the eigenvalues of the matrix 

R(t)=p(t)p(t), (43) 

with p(t) given by pit) = a y £g> o~ y p*it)o y £g> o~ y , o y being the Pauli matrix, and pit) the 
density matrix representing the quantum state of the system. Inserting Eq. fj4Tl) into 
Eq. ( j4"3j) . it is possible to derive the time evolution of concurrence, studying in this 
way the features of the non-classical correlations characterizing the system of the two 
coupled qubits. More in details, it is possible to see that when the system is initially 
prepared in the factorized state 1 01) its dynamics is characterized by the existence of 
damped oscillations corresponding to the periodic appearance and disappearance of non- 
classical correlations between the two qubits. The oscillations in the concurrence are a 
direct signature of the well known oscillations of the single excitation between the two 
qubits, which are caused by the resonant terms in the qubit-qubit interaction. On the 
other hand the stationary entanglement at long times is a consequence of the fact that, 
due to the counter- rotating terms in the Hamiltonian in Eq. (TjQ), the ground state of 
the system is not the unexcited state 1 00) , but the entangled state \a) in Eq. ([3]). From 
Eq. (00) it is straightforward to see that the smaller the coupling constant A the smaller 
the amount of stationary entanglement. 

The system dynamics is even richer when the system is initially prepared in a state 
with two excitations or in an arbitrary superposition of the states 1 00) and |11), i.e. 



|V>(0)> = v/p|00> +e l Vl-p|ll>. (44) 
Also in this case p is a non-negative real number < 1, while we have: 



Paa(0) 
Pdd(0) 



VP cos (^j - e^y/l-p sin (^j 



VP sin ( ~f) + e^V 71 ~P cos (tJ 



(45) 



(46) 



PaM = L/P cos r-± \ - e^VW sm f | 



x ( VP sin ( |] + e-^V^P cos (| ) ) > (47) 
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with pa = and pij = otherwise. With such an initial condition, the density matrix 
of the system at a generic instant of time assumes the form 

p(t) = Paa(t)\a)(a\ + p bb (t)\b)(b\ + Pcc(t)\c)(c\ + Pdd{t)\d)(d\ 

+ p ad (t)\b)(c\+p da (t)\c)(b\. (48) 

Exploiting Eq. (jUj) (with p = 0) as well as Eqs. (14"5]) - (14"8"]1 it easy to derive the 
time evolution of the concurrence of the two qubits. Figure 3 shows that, apart from 
the oscillations of the entanglement due to the oscillations of the excitation in the one- 
excitation subspace, the behaviour of the concurrence is characterized by the phenomena 
of sudden birth and sudden death of the entanglement already found by Ficek and Tanas 
P3] for the scenario of two interacting qubits. While the phenomenon entanglement 
sudden death, i.e., the complete disentanglement of the system in a finite time, is 
quite well understood and is a feature common to any dissipative two-qubits dynamics, 
provided the system starts from the proper initial state [TTJ [12], the occurrence of 
stationary entanglement is a phenomenon usually ascribed to the presence of a reservoir 
which is common to the two parts of the bipartite system [TH [16]. What we have 
proved here is that the same phenomenon can occur even in the presence of independent 
reservoirs for the two qubits, provided the interaction Hamiltonian between the two 
qubits contains also the counter-rotating terms which are usually neglected in the study 
of the dynamics. This is due to the fact that, in the weak damping limit, the quantum 
jumps occur between eigenstates of the system Hamiltonian and to the fact that, due 
to the presence of counter-rotating terms in the two-qubit interaction, the ground state 
of the system is an entangled state, as we have seen before. 




t 

Figure 2. Concurrence of the two qubits as a function of time t (in units of 10*(A) _1 ). 
The initial state is |01), the parameters assume the values wi/A = LOij\ = 10, 
Ti = T 2 = K, 7j,h/A = 7//41/A = 7/,22/A = 7/,22/A = 0.01. 
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5. Discussion and Conclusive Remarks 

To summarize, we have presented a microscopic derivation of the Markovian master 
equation in the weak damping limit governing the time evolution of two interacting 
two-state systems, each of them coupled to independent bosonic reservoirs, and we have 
studied the time evolution of the entanglement between the two qubits for various initial 
conditions of the bipartite system. 

The behaviour of the concurrence is consistent with the results of the literature, 
showing the occurrence of entanglement sudden death for some initial states and of 
stationary entanglement for any initial state. The presence of stationary entanglement 
which is not washed away from dissipation is perhaps the most important result of our 
analysis. Indeed this feature is usually ascribed to the presence of a common reservoir 
which correlates the two parts of the bipartite system. Instead we have shown here that 
stationary entanglement can occur also due to the interaction between the two parts, 
provided their interaction is described in a complete way which includes the counter- 
rotating terms in their interaction Hamiltonian. We feel that our approach could help 
in understanding the role of energy non-conserving terms in the dissipative dynamics 
of open bipartite quantum systems. From this point of view more is expected from the 
analysis of the cases with reservoirs at different temperatures or with non-flat spectrum, 
in the framework of a non-Markovian extension of the master equation presented here. 
These points will be the subject of our future work. 
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Figure 3. Concurrence of the two qubits as a function of time t (in units of 10*(A) 1 ). 
The initial state is |11), the parameters assume the values u>i/X = uj 2 /X = 10, 
Ti = T 2 = K, 7/,11/A = 777,n/A = 7/,22/A = 7/,22/A = 0.01. 
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